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Abstract 

The viscoelastic material functions for the Becker and the Lomnitz 
rheological models, sometimes employed to describe the transient flow 
of rocks, are studied and compared. Their creep functions, which 
are known in a closed form, share a similar time dependence and 
asymptotic behavior. This is also found for the relaxation functions, 
obtained by solving numerically a Volterra equation of the second kind. 
We show that the two rheologies constitute a clear example of broadly 
similar creep and relaxation patterns associated with neatly distinct 
retardation spectra, for which analytical expressions are available. 

1 Introduction 

The purpose of this paper is to draw the attention of polymer scientists 
on two models used in Earth rheology. They are usually refereed to as 
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Becker and Lomnitz to honor the scientists who have introduced them in 
1925 [T] and in 1956 [2], respectively. Both models exhibit slow varying creep 
laws suitable for simulating the flow and the (quasi frequency independent) 
energy dissipation in rocks, see e.g. Strick and Mainardi |3]. Though the 
corresponding relaxation laws are not considered in geophysical frameworks, 
they are certainly of interest in the theory and applications of linear 
viscoelasticty. As far as we know, in both classical and contemporary polymer 
science, see e.g. [H El [6l [7], these rheological models have not been taken 
into account. 

In the following we will discuss the analytical creep laws for the two 
models along with their graphical representation versus dimensionless time 
both in linear and logarithmic scales. Because the differences between the 
two creep laws remain small as time is evolving, we also show the rate of creep 
in order to have a better insight of the comparison. Then, we numerically 
compute and visualize the corresponding relaxation laws by solving a Volterra 
integral equation of the second kind. The major difference between the two 
models is found in their retardation spectra. 

2 The creep laws 

In Earth rheology, the law of creep is usually written as 

J{t) = Ju[l + qm]. t>0, (1) 

where t is time, Ju is the un-relaxed compliance, g is a positive dimensionless 
material constant, and i/j{t) is the dimensionless creep function. Consistently 
with the general theory of linear viscoelasticity, ipit) is a Bernstein function, 
that is positive with a completely monotone derivative, with a related 
spectrum of retardation times (see e.g. Mainardi [8]). 
For the Becker model pLj we have 

7/'^(t) = Ein(t/ro) , t>0, ro>0, (2) 

where Ein denotes the modified exponential integral function (see e.g. [8]). 
Assuming tq = 1, we have the integral and series representations 

/•* 1 _ p-U ^ +n 

Em(t)= du = Y(-ir-^ -, t>0, (3) 

Jn u ^ nn\ ~ ^ ^ 

" n=l 



2 



t Log(t) 

Figure 1: The creep functions for the Becker (B) and Lomnitz (L) rheologies, 
shown clS 8b function of time. 

hence the rate of creep is 

dt ^ ' t Mn + 1 ! ' - ^ ^ 

71=0 ^ ' 

For the Lomnitz model [2] we have 

V^^(t) =log(l + t/ro), t>0, ro>0, (5) 

where log denotes the natural logarithmic function. Taking again tq = 1, we 
have the series representation 

iog(i+t) = V(-ir-^-, t>o, (6) 

^-^ n 

n=l 

which implies a rate of creep 

n=0 



For the Lomnitz law the series representations given by Eqs. (6) and (7) 
are convergent only for < t < 1, at variance with Eqs. (3) and (4) for 
the Becker law that are convergent for all t > 0. However, all these power 
series are suitable only for sufficiently small times because their numerical 
convergence falls down very soon. 
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Figure 2: Rate of creep for the two models, as a function of time. 

In order to compare the creep behavior of the two models, we show 
[t) , ip^ (t) and their time derivatives in Figures 1 and 2, respectively, taking 
both linear frame (a) and logarithmic (b) time axes. The overall similarity 
between the two creep functions is apparent in Figure 1, also showing that 
the Becker rheology accounts for a somewhat larger strain relative to Lomnitz 
(i.e., > ^^)- Inspection of the corresponding rates of creep in Figure 2, 
clearly shows that, for finite values of time, the Becker creep systematically 
evolves at a larger rate with respect to Lomnitz (i.e., dip^ /dt > dip^/dt). For 
long times, both rates of creep decay to zero as 1/t as we easily note from 
Eqs. (4) and (7). 

3 The relaxation laws 

The relaxation modulus G{t) for the two rheological models can be derived 
from the corresponding creep laws through the general Volterra integral 
equation of the second kind [8] 



As a consequence, the dimensionless relaxation function defined by 





cj){t) = JuG{t) 



(9) 



obeys the integral equation 




(10) 
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where the rate of creep is given by Eqs. (4) and (7) for the Becker and the 
Lomnitz laws, respectively. In order to solve numerically Eq. (10), we have 
used standard numerical methods. 
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Figure 3: Relaxation functions for the two rheological bodies, as a function 
of time. 

The results are shown in Figure 3, assuming q = 1 and adopting 
again both linear and logarithmic time axes. As expected according to the 
similarity of the corresponding creep functions, the two relaxation functions 
show similar features for the two models. However, it is apparent that 
the Becker model exhibits, at a given time, a somewhat larger amount of 
relaxation relative to Lomnitz (i.e., (p^ < (p^). By visual inspection of the 
curves it also appears that the Becker rate of relaxation exceeds that of 
Lomnitz model. 



4 The retardation spectra 

The determination of the time-spectral functions from the knowledge of the 
time-dependent material functions J(t) and G{t) is a fundamental problem 
from theoretical and experimental view points in polymer science. It can be 
formally solved through an analytical method outlined by Gross |9] based on 
the Laplace transform of the time derivatives of J(t) and G(t), see also [8]. 
For the present models we use the Gross method to derive the retardation 
spectrum from the Laplace transform of the rate of creep. By definition of 
the retardation spectrum Re{T), we have 

J{t) = Ju + j^ Re{r) (l - e"^/^) dT , (11) 
where r denotes the retardation time. 
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For the two rheological models considered in this note, closed-forms exist 
for the retardation spectra. For the Becker model, according to [9], the 
retardation spectrum turns out to be discontinuous, with 



R?(t) = -H(t-1] 



[12) 



where H{t) denotes the Heaviside step function, while from \T0] for the 
Lomnitz model we have the continuous spectrum 



T 



(13) 



The two spectra are compared in Fig. 4 as function of r. Although these 
spectra show a dramatic difference character, we note that they both show 
a peak for r = 1 and that they both decay, for r oo, as 1/r. 
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Figure 4: Retardation spectra for the two rheological models as a function 
of the retardation time. 



5 Conclusions 

In this paper we have discussed and compared the time-dependent material 
functions for two viscoelastic models introduced by Becker and Lomnitz 
known in Earth rheology. These functions for the two models show a 
broadly similar behavior, so that they could hardly be discriminated from 
experimental point of view. Despite this similarity, the corresponding 
retardation spectra show a dramatic difference. While the Lomnitz spectrum 
varies smoothly on the whole range of retardation times, the Becker one 



6 



displays a cut off for short time even allowing a similar decay at large 
times. The examples discussed in this paper clearly show that in order to 
discriminate between two rheological models exhibiting very similar creep 
and relaxation behaviors the evaluation of the spectra is required. 
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